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^ \ We present a novel theory of charge-carrier mobihties in organic molec- 

O 

ular crystals of high purity. Our approach is based on Holstein's original 
<^ ^ concept of small-polaron bands but generalized with respect to the inclusion 

' of nonlocal electron-phonon coupling. We derive an explicit expression for 

in ' 

\ the mobilities as a function of temperature and, using ab-initio methods to 

I obtain the material parameters, we demonstrate its predictive power by ap- 

\ plying it to naphthalene. The results show a remarkably good agreement with 

, experiments and provide new insight into the difference between electron and 

Q I hole mobilities as well as their peculiar algebraic and anisotropic temperature 

O ' 

dependences. 
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Typeset using REVT^ 



Organic semiconductors with 7r-conjugation are very promising materials for low-cost 
and easy-to-process electronic and optoelectronic devices such as light-emitting diodes [1,2], 
lasers [3,4], solar cells [5,6], and thin- film transistors [7-11]. Besides the polymeric materi- 
als, a very important class of organic semiconductors are molecular crystals of high purity. 
Due to their high degree of structural order, such crystals are ideal candidates for the in- 
vestigation of the intrinsic excitations and charge-carrier transport phenomena in organic 
solids. Consequently, many interesting fundamental experimental and theoretical studies 
have been performed in recent years and new results emerge at a rapid pace [12-16]. A 
particularly important and challenging topic is the understanding of the charge- carrier mo- 
bilities in these crystals, as the mobility is a fundamental material property and a central 
quantity for optimization of device performance. 

In comparison to covalently bonded inorganic semiconductors, organic molecular crys- 
tals exhibit weak intermolecular bonds and, hence, narrower electron bands and stronger 
electron-lattice interaction. As a result, the dressing of the charge carriers by phonon clouds, 
i.e., polaron formation, becomes quite important. Pioneering theoretical work in this area 
was done by Holstein [17] who introduced the concept of small-polaron bands for a local 
(on-site) electron-phonon coupling. The main prediction of his theory is that the mobilities 
initially decrease with rising temperature T due to bandwidth narrowing but increase again 
for higher T due to phonon-assisted hopping. Two decades after Holstein's prediction, the 
interplay between metallic (bandlike) conduction and activated (hopping) transport in or- 
ganic crystals was reported for the electrons in naphthalene but not for the holes [18-20]. 
Later, Kenkre et al. were able to fit the measured electron mobilities reasonably well to 
Holstein's model, assuming directionally-dependent local-couphng constants [21]. Despite 
the success of such a fitting procedure, the lack of a first-principles description of charge- 
carrier mobilities in organic crystals has left several fundamental questions unanswered. 
This concerns, in particular, the different behavior of electrons and holes, the microscopic 
origin of the crystallographic anisotropy in the T dependence, and the infiuence of nonlocal 
(Peierls-type) couphngs such as present in the Su-Schrieffer-Heeger model [22]. Moreover, 
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the fundamental question has been raised whether or not the concept of small-polaron bands 
is apphcable at all, since band motion may be washed out by even weak disorder [23]. 

In this Letter, we answer the above questions by developing a novel first-principles theory 
of charge-carrier mobilities in organic crystals. We present strong evidence for the usefulness 
of the small-polaron band picture once Holstein's original model is extended towards the 
inclusion of nonlocal electron-phonon coupling. While a few theoretical studies on nonlocal 
coupling have been performed in the past [24-28], it was only recently that the importance 
of these contributions could be demonstrated explicitly in the case of polaron bandwidth 
narrowing in oligo-acene crystals [29]. Here, by means of the Kubo formalism, we obtain 
an explicit expression for the mobilities, including both local and nonlocal coupling. The 
practical usefulness of this approach is then exemplified by ab-initio studies for naphthalene 
(CioHg) crystals. 

We consider a mixed Holstein-Peierls model for the interaction between electrons (holes) 
and phonons. In a tight-binding description, this corresponds to a Hamiltonian of the form 

H = Hei + + -ffei-ph where 

mn 

i^ph = E^^Q(^Q^Q+i) ' (2) 
Q 

-ffel-ph = Yl ^'^Q5'QmnaL«n(&Q+&-Q) , (3) 
mnQ 

describe the electrons, the phonons, and the electron-phonon interaction, respectively. Here, 
the operators and b^Q := 6^^ annihilate (create) an electron at equivalent lattice sites 
= {Ri, . . . , Riv} with energy and a phonon in the mode A with wavevector q and 
frequency cux^ —: cuq, respectively. The strength of the electron- lattice interaction in Eq. (3) 
is determined by the dimensionless quantities gQmn, which describe both the /oca/ coupling to 
the on-site energies Emm (Holstein model) and the nonlocal coupling to the transfer integrals 
£mn {iTi 7^ n, Pcicrls model) . The latter type of coupling will be shown to play a crucial role 
in the understanding of the anisotropic T dependence of the mobilities. 

For the description of the electrical conductivity, we use linear-response theory and intro- 
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duce a conductivity tensor aa/siuj) that relates the electric field E(a;) present in the sample 
to the current 3{u)) that is induced by this field, Jq,(u;) = J2i3<^a/3{^)Ei3{u). Here, we con- 
centrate on the mobilities /x^ which are obtained in the zero-frequency limit according to 
Caa(<^^0) = NcCo/ia, whcrc cquals the negative (positive) elementary charge in case of 
electron (hole) conduction and Nc denotes the number of charge carriers. In thermal equi- 
librium, the mobilities Ha can be calculated by means of the Kubo formula for electrical 
conductivity [30], 



which relates the mobility to a microscopic current-current correlation function. The r.h.s. 
of Eq. (4) contains the quantum-mechanical current operator j{t) — e^^* j e"^^* where j is 
obtained from the polarization operator P = cq Y,m R-mflm'^m Heisenberg equation 

of motion j = ^ = ^^[P, It is important to stress that the current j consists of two 
different contributions, j := j^''^^ + which are given by 



operators, the second term j^^^^ describes a phonon-assisted current that stems solely from 
the nonlocal electron-phonon coupling and has no counterpart in local-coupling theories. 

The evaluation of the Kubo formula (4) is performed by the method of canonical transfor- 
mation [30] within a generalized version that accounts for both local and nonlocal coupling 
[27-29]. By means of this technique, the electron-phonon interaction is incorporated non- 
perturbatively, in accordance with the underlying physical picture of small-polaron bands. 
Here, we skip the details of the quite tedious derivation and the approximations involved 
(which will be pubhshed elsewhere [31]) but rather concentrate on the final results for the po- 
laron mobilities in the case of dispersionless optical phonons (cuq — > 0;^) for which symmetry 
requirements impose the form gQrnn = ^:^5'Amn(e"''^'^'"+e"*'i'^"). 




(4) 



mnQ 

and originate from ifei and i/ci-ph, respectively. While the current j*^^^ contains only electron 
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It is instructive to consider first the mobility /i^J^ due to the currents j^^^ in Eq. (5), 
which we obtain as 

2kBin 7-00 

where A^a= (e^^/*^^-^— denote the phonon occupation numbers. The extra factor involv- 
ing the phenomenological broadening parameter F accounts for scattering processes beyond 
the model, e.g., due to defects. Those terms in the exponent of Eq. (7) that contain the time- 
dependent auxiliary function ^\{t) = {1+Nx)e''"^^^ + Nxe^'^''^ describe mco/ierent scattering 
events involving actual changes in phonon numbers (hopping) whereas the remaining terms 
account for purely coherent scattering processes (bandwidth narrowing) . While our above 
result for fi^'^ can be compared to the Holstein-model mobilities obtained by Kenkre et al. 
[21,32], an important difference is that the exponent in our formula (7) is not only governed 
by the local coupling but by effective coupling constants G\ — {gxmmY + \ Hki^miQxmk)'^ 
that are composed of both the local and nonlocal ones. 

So far, we have neglected the contributions from the phonon-assisted currents j^^^-' of Eq. 
(6). If we include them, we obtain the total mobilities /Xq within the Holstein-Peierls model, 
which can be expressed in a form identical to Eq. (7) for but with the replacement 

{Smnf {emn-^mnf + \ Y.i^^ ^9 Xmnf ^ x{t) ■ (8) 

X 

Here, we have introduced the 

shorthand notation A^„=| Ea T^^x [9xmn{9xmm+9xnn) + \ T.k^m,n 9xmk9xkn\ ■ An essential 
qualitative difference between ji^^^ and is that only the total mobilities //^ can account 
for an anisotropic T dependence. This is due to the additional terms containing ^\{t) in Eq. 
(8), which describe incoherent scattering processes that are absent in local-coupling theories 
and reflect an inherently transport-promoting effect solely caused by the nonlocal coupling. 
The strength of the resulting phonon-assisted hopping and, hence, the anisotropy in the T 
dependence is basically determined by the ratios huJxgxmn/^mn, which may be dramatically 
different for different directions. 



As a first appfication of the above theory, we perforin model studies for naphthalene, 
which crystallizes in a monoclinic structure (P2i/a) and exhibits a herringbone stacking with 
two equivalent molecules per unit cell. In order to obtain the material-specific parameters 
R-m, £mn, gxmn, ^ud cux we use the following 3-step strategy. 

First, we determine the equilibrium structure of the crystal by means of state-of-the-art 
DFT-LDA calculations using the ab initio total-energy and molecular dynamics program 
VASP [33]. For the lattice parameters we obtain the values a — 7.68 A, b — 5.76 A, c = 
8.35 A, and the monoclinic angle P — 125.7°. For this geometry, the mtermolecular optical- 
phonon energies {hcui = 10.7 meV, hu!2 = 14.2 meV, hcu^ = 17.4 me V) and polarizations 
Bx are obtained within the rigid-molecule approximation, using a doubled Brillouin zone 
corresponding to the lattice {Rm} of all the molecules [34]. In principle, our approach 
does also allow the inclusion of intramolecular phonons. However, in a recent Letter on 
oligo-acene molecules [12] it was reported that those intramolecular phonons that couple 
most strongly to the electrons have significantly higher frequencies than the intermolecular 
phonons considered here. Hence, while the zero-point fluctuations of the intramolecular 
phonons will lead to an isotropic overall reduction of the mobility values [35] , their influence 
on the actual T dependence will be small due to the low occupation of these modes up to 
room temperature. 

In the second step, we obtain the values from a flt of the ground-state ab ini- 
tio HOMO and LUMO energy bands to a tight-binding model, including the on-site 
energy and the six most important transfer integrals between nearest neighbors, i.e., 
{■mn}—{0, a, b, c, ac, ab, abc} belonging to R^— R„ = 0, ±a, ±6, ±c, ±(a-|-c), ±(|±|), and 
± ( I ± I + c) , respectively. 

In the third step, we rotate the molecules by amplitudes Aux according to the polar- 
izations ex of phonon mode A ("frozen phonon") and flt the new resulting ab initio band- 
structure again to the tight-binding model. Then, the electron-phonon coupling constants 
gxmn are obtained from the changes in the transfer integrals by numerical differentiation 
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ffAmn^fi— ^^^" ; in the limit Aux — > 0. All calculated parameters are compiled in Table 
I. Note that our calculated effective coupling constant (5'cfr)^:= Z]a=i ^a=(1-52)^ for the 
LUMO band is comparable in value to the directionally-dependent fitted coupling values of 
1.62, 1.83, and 1.88, as obtained by Kenkre et al. [21]. Finally, for the line broadening, we 
have chosen a small value of = 0.1 meV, corresponding to the case of ultrapure crystals 
[36]. 

In Fig. 1, we present the electron and hole mobilities in naphthalene crystals as a function 
of temperature T for the a, b, and c' directions, with c' being perpendicular to the a6-plane 
of the molecular layers. We make three important observations. First, the hole mobilities 
are generally larger than the corresponding electron mobilities. Only in the a direction do 
we observe comparable mobilities for electrons and holes at T fa 300 K. This general trend is 
in good agreement with experiments of Karl and co-workers [19,20]. Importantly, the lower 
electron mobihties cannot be understood by simply looking at the transfer integrals of Table 
I but they are mainly due to their stronger couphng to the phonons (LUMO: Qes — 1.52, 
HOMO: g^e = 0.92). 

Second, while the hole mobilities in the different directions decrease very similarly with 
T, the electron mobilities exhibit a pronounced anisotropic T dependence. In particular, we 
find a metallic (bandlike) behavior within the a6-plane and a slightly activated (hopping) 
transport in the c' direction. Again, these findings from our ab initio calculations are in 
good agreement with experimental data [18,19], which show such a nearly T-independent 
mobility around T pa 150 K exclusively in the c' direction. As discussed above, this pro- 
nounced anisotropy is due to the phonon- assisted currents j*^^^^ and depends on the ratios 
Ty^\g\mn/ ^mn- For the electrons, this ratio is particularly high for hcu2g2c/^c, which explains 
the peculiar flattening of /^c' that is not seen at all for /x^P alone. Thus, our findings show 
clearly that anisotropy effects in the T dependence can be explained in a natural way by 
the inclusion of nonlocal electron-phonon coupling, and do not require the use of additional 
concepts such as a directionally-dependent local couplings [21] or correlated successive hops 
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of strongly localized charge carriers [37] . 

Third, our calculated hole mobilities obey at elevated temperatures the power law 
Ha oc T~^-^. This agrees nicely with Karl's experiments [19] who found iia oc T~'^ (7 = 2.9, 
2.5, 2.8). This peculiar power-law dependence, which has been poorly understood so far and 
for which several explanations have been put forward [20], can be straightforwardly under- 
stood from our theory. In the limit of small broadening hV and for /Xq ~ A*^^-* (here fulfilled 
for the holes), a spectral analysis of Eq. (7) shows that the dominant electron- phonon scat- 
tering processes are those that allow only energy exchange within each phonon mode but not 
between different modes, and the mobility can be expressed in terms of modified Bessel func- 
tions, //i^) oc (nLi e-2«>(i+2iv,)j^(^^)j^ ^^ere zx = AG^^N^il+Nx). The limit of high 
T is obtained using Io{zx^oo) ^ and we find //(f) oc T'^ (ni=i ^"°-^) = T'^'^ 

which indicates that the exponent is actually a measure for the number of relevant phonon 
modes per molecule. For the electrons, the situation is more complicated due to a subtle 
interplay between the now significant contributions of j^^^^ and the larger g^s value. This 
prohibits an easy analytical evaluation but, by looking at our calculated in-plane electron 
mobilities in Fig. 1, we find at high T an approximate T~^-^ dependence, again in close 
agreement with the experimental observations [19,20]. 

In summary, we have presented a novel theoretical description of charge- carrier mobilities 
in ultrapure organic crystals based on a mixed Holstein-Peierls model where not only the lo- 
cal but, importantly, also the nonlocal electron-phonon coupling is included. The predictive 
power of our theory and the usefulness of the small-polaron band concept have been demon- 
strated by ab initio studies for naphthalene crystals where our results show a remarkably 
good agreement with experiments and provide new microscopic insight into several hitherto 
poorly understood phenomena including the different behavior of electrons and holes as 
well as the peculiar algebraic and anisotropic temperature dependences. Finally, since our 
theory may also be applied to other weakly-bonded materials with strong electron-phonon 
coupling, it may become an important contribution to the understanding of the intrinsic 
charge-carrier transport in these materials, too. 
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TABLE I. Top: 


transfer 


integrals Smn 


and electron-phonon 


coupling 


constants 


9Xmn for 



the HOMO and LUMO bands of naphthalene crystals. Bottom: effective coupling constants 

(ffefr)':=EA=i Ga, where Ga = 9lo+9la+9lb+9lc+9lac+^9lab+^9labc- 
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FIGURES 

FIG. 1. Charge-carrier mobilities vs temperature T for the a, b, and directions in naphthalene 
crystals. Top 3 curves: /Uq, for holes, middle 3 curves: for electrons, bottom curve: for 
electrons. As a guide to the eye, several power laws /x^ oc T~"' are also depicted. All mobilities are 
plotted in units of x ^^^^ but an unknown prefactor r [35] should be included. 
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